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A description of low-cncrgy scattering of pions and nuclei within a BUU transport model is 
presented. Implementing different scenarios of medium modifications, the mean free path of pions 
in nuclear matter at low momenta and pion absorption reactions on nuclei have been studied and 
compared to data and to results obtained via quantum mechanical scattering theory. We show 
that even in a regime of a long pionic wave length the semi-classical transport model is still a 
reliable framework for pion kinetic energies greater than 20 — 30 MeV. Results are presented on 
vr-absorption cross sections in the regime of 10 MeV < T^^„ < 130 MeV and on photon-induced 
TTTT production at incident beam energies of 400 — 500 MeV. 
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I. INTRODUCTION AND MOTIVATION 

According to our present knowledge, the properties of hadrons and their interactions can be described in terms of 
Quantum Chromo Dynamics (QCD). In the limit of vanishing quark masses QCD incorporates chiral symmetry, which 
is spontaneously broken in vacuum. The order parameter (qq) of this symmetry breaking is expected to decrease by 
about 30% already at normal nuclear matter density Therefore, signals for partial chiral symmetry restoration 

should be observable in nuclear reaction experiments and, in particular, photon induced processes are highly suited 
due to two key reasons. The reaction leaves the nucleus close to its ground state, so the reaction takes place under 
well defined conditions. And as a second point, the photon penetrates deeply into the nucleus, giving rise to a high 
effective density. 

The modification of the so-called a or /q (600) meson inside the nuclear medium was proposed as a signal for such a 
partial symmetry restoration. Theoretical models predict a shift of its spectral strength to lower masses and a more 
narrow width due to the onset of the restoration (0; 0). The cr- meson is a very short-lived state with a width of 
roughly 600 — 1000 MeV (@), decaying predominantly into a tttt final state in S'-wave. Therefore, the experimental 
aim has been to find modifications of this state in tttt production reactions in finite nuclear systems close to threshold. 

Such experiments have been performed with incident pions by the CHAOS collaboration (0; @) and with photons 
by the TAPS collaboration |2; Jil; |^; IT^ . Both experiments have shown an accumulation of strength near the tttt 
threshold in the decay channel of the a in large nuclei. A possible interpretation of this effect is the in-medium 
modification of this resonance due to partial symmetry restoration. Already in an earlier work ijlST) . we have, however, 
pointed out the importance of conventional final-state effects in the analysis of the experiment performed by the TAPS 
collaboration. 

To analyse the experiment performed by the TAPS collaboration (0; 0; ^2 0)i we need to consider pions with 
very low energy in nuclear matter. Aiming to employ a BUU transport model, we first need to establish that it is 
meaningful to discuss pions with a long wave length in such a semi-classical picture. Actually, this question can only 
be answered by comparison to quantum calculations and experimental data. Therefore, in this paper we investigate 
the validity of our model for the treatment of the tttt final state with the description of pion-induced scattering and 
absorption experiments within the BUU-framework. 

Alr eady in earlier works of Salcedo et al. (1141) . with a simulation of pion propagation in nuclear matter, and of Engel 
et al. lfl5() . with a precursor of our present simulation, pions with kinetic energies in the regime of 85— 300 MeV have 
been investigated in transport models. As motivated above, we now investigate even less energetic pions. Therefore, 
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we have carefully accounted for Coulomb corrections in the initial channel of 7r-induced processes and improved the 
description of the threshold behavior of the cross sections in the model. 

This article is structured in the following way. First we introduce our Giessen BUU (GiBUU) transport model 
for the treatment of photon- and pion-induced nuclear reactions emphasizing the most relevant issues. Next, we 
present results on the pion mean-frce-path in nuclear matter and discuss the consequences of medium modifications 
of the pion. Hereafter we compare our simulation to experimental results on pion scattering off complex nuclei to 
validate our model. Finally, we conclude with results concerning the tttt photoproduction at incident beam energies 
of 400 - 500 MeV. 



II. THE GIBUU TRANSPORT MODEL 

Boltzmann-Uehling-Uhlenbeck (BUU) transport models are based on the Boltzmann equation, which was modified 
in the 1930's by Nordheim, Uehling and Uhlenbeck to incorporate quantum statistics. It is of semi-classical nature, but 
can be derived as a gradient-expansion of the fully quantum Kadanoff-Baym-equations |)16() . Nowadays this model is 
widely used in heavy- ion collisions, electro- and pion-induced reactions to describe these processes in a coupled-channel 
approach. 

The Giessen BUU-model has been utilized in various applications including heavy-ion-collisions, photon-, electron-, 
pion- and neutrino-induced processes. In contrast to earlier works (ll.^ i llTtlTsI) we are now using a new implementation, 
named GiBUU ITol) . It is in its algorithms based upon the old implementation, therefore also reproduces the former 
results. For a detailed discussion concerning the physical input and the algorithms we refer the reader to ifl^ IT8 H20|) 
and the references therein. For the reader interested in technical details, we note that the new code is now written 
in modular manner in modern Fortran l)2l|) utilizing Subversion l)22() version control in a multi-user environment. It 
represents a major leap forward to a simplified handling of the transport code and a more transparent implementation 
of future projects. 



A. The BUU equation 

The BUU equation actually consists of a series of coupled differential equations, which describe the time evolution 
of the one-particle phase-space-densities /f (p, r, t) in the limit of low particle densities. The index a = t:,uj, N, A, . . . 
denotes the different particle species in our model. 

For a particle of species X with negligible width, its time evolution is given by 

d/f (f , p, t) _ Off (f, p, t) , dHx a/f (r , p, t) OHx df^ (f , p, t) _ ^ (rx ra rb ^ r,. 
dt - dt ^ dp dr dr dp -^-oii[Ji ,h,h,---) 

with the one-body Hamilton function 



Hx{r,p) 



\j(p + Ax (7% p))" +m\ + Ux (r, p) + A\ (f, p) 



The scalar potential Ux and the vector potential A^ of species X may in principle depend upon the phase space 
densities of all other species. Hence, the differential equations are already coupled through the mean fields. In the limit 
of IcoU = 0, equation ^ becomes the well-known Vlasov equation. The collision term Icon ({/f , a = n,uj, N, A, . . .}) 
on the right-hand side of equation ^ incorporates explicitly all scattering processes among the particles. The reaction 
probabilities used in this collision term are chosen to match the elementary collisions among the particles in vacuum. 
Within the BUU framework the total reaction cross sections are given by an incoherent sum over all partial cross 
sections. Interference effects are therefore neglected. 



B. Elementary processes 

The elementary processes of interest in the present study are ttN — * ttN and ttNN ^ NN reactions. Both 
are modeled within a resonance picture. Therefore, we express the relevant cross sections as an incoherent sum of 
resonance and background contributions. 
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The resonance cross sections are obtained from the partial wave analysis of ref. lj23|) . The background cross sections 
denoted by a^^ are chosen in such a manner, that the elementary cross section data in the vacuum are reproduced (|20|) . 
Within our model, background contributions are instantaneous in space-time, whereas reactions through resonance 
channels are delayed by the lifetime of the resonance. 

All different events in pion nucleon scattering can be categorized into four different isospin channels 



^ TT n — >7r n '-'"tt+p— »7r+p ; 



^tt'-^P^tt^P • 



(3) 
(4) 
(5) 
(6) 



The cross sections in the individual channels arc cither connected by time reversal or isospin symmetry. The first 
channel (eq. is a pure isospin 7 = 3/2 scattering process, whereas the other three channels are mixtures of / = 1/2 
and / = 3/2. The cross section for the / = 3/2 channel CTtt n—^a^tt n is given explicitly in |23) based on the resonance 
analysis by Manley and Saleski |23l) . As can be seen in the upper left panel of fig.H the first channel is well described 
solely by the A-resonance scattering; therefore, we neglect a background term in this channel. In the second and third 
channels (eq. ^andEJ there are good data sets for 



TT p 



• TT p 

TT^ n. 



down to very low energies. Hence, we introduce a background term on top of our resonance contributions for a better 
description of those channels as shown in the remaining panels of fig. The last channel (eq. is inaccessible for 
experiment, therefore we can not introduce any background term. 

Concerning total cross sections, the relevant contributions to the scattering of pions and protons are shown in fig. El 
for the different pion charge states. The cross sections on the neutron follow by isospin inversion. 
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FIG. 1 Cross sections for elastic scattering and charge exchange reactions in pion-nucleon scattering. The solid lines show the 
full cross sections, whereas the dashed lines represent the background contributions. Data are taken from lj2!Tl : \2(t\ . 
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FIG. 2 Elementary cross sections for elastic and inelastic scattering of pions and protons. Note that the resonance production 
processes contribute to elastic and inelastic scattering processes. The data ||2RI: l2^: l23l are for a"(tot). 



The process NN NNtt has been extensively studied in several experiments over the last twenty years (|26l : 

S E3 Sil El S El 

We can therefore construct well defined background cross sections on top of our resonance 
contributions for all possible isospin channels. In fig.|2|we show the relevant cross sections. 

The ttNN — > NN process is described by the two-step process wN R followed by RN NN, and a direct 
background contribution. One defines the pion absorption rate 

T. c c IPabI \M\'^ PNa PNb (r,. 

r^,Ar,._^„Ar, - . (7) 

This rate depends on the densities pma ^-iid pnb of the nucleons in the initial state. The matrix element Ai can be 
calculated by detailed balance, assuming that it only depends on the Mandclstam variable s. 

|A^(s)| = \MN^Ni^NANBTr{s)\ \MNANBTr~^ N^N^s)] 

To obtain this matrix element we consider more closely the NN NNtt process. The cross section for this process 
is given by 

, I ^ I — -^S"^ {pa + Pb - PA — PB — Ptt) IMNaNb-'NANBTr]'^ (8) 

d^PA d^PB d^p„ 



(27r)3 2Ea {2tt)32Eb {2T:f2E^ 



1 /•(miB)r,.ax /■(mi„)max 

Sab———--,\M\' / dmls / dmL (9) 



(27r)364|p;b|\/s J(^n\^)n.^^ Arn 
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FIG. 3 Elementary cross sections for diflterent A'^ — * NNn isospin channels. The solid lines show the full cross section, whereas 
the dashed lines represent the non-resonant background contribution. The data are taken from ||2^ : l29t Isot Isil |33 133: |3^ . 
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as symmetry factor; pat denotes the momentum of the particles a and b in their center of mass frame. Equation (PJ 
holds under the assumption that the matrix element is only dependent on the Mandelstam s. So we get 

\M{s)\^ = isAB———^———3 I dm^gdm^A aN^N.^N^NB-^ ■ 
\ 64(27r)'^ \pab\ \/s J J 

After inserting this result for |A^(s)p into equation {Tj) we find that rjv^Ar^Tr-tAfaJVt, depends linearly on (JNaNb^NANB-n- 
This cross section CFNaNb^NANsTr is, according to equation ((21), a sum of background and resonance contributions. 
Therefore, ^NANsT^^NaNt can also be split into a resonance and a background contribution. 

T^BG I prcsonancc contribution 

J- NANgTr^NaNb — -L NANsTr-^N^Nh ' NANBTr-*N^Nb 
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FIG. 4 The NNn NN background absorption rate for symmetric nuclear matter at p = po . In-medium-modifications are 
not included in these decay widths. The rate for n" is shown on the left and for vr on the right panel. 



with 

presonance contribution resonance contribution 

More details concerning this absorption rate can be found in ijlTfl and l)20l) . The resonance absorption part is 
included in the collision term for the resonances, which are propagated explicitly. The background absorption rate 
^WANsTr—fNaNb shown in fig. 01for symmetric nuclear matter at p = po with po ~ 0.168 fm~"^ being normal nuclear 
matter density. For positive pions we get the same results as for negative ones due to isospin symmetry. Notice that 
in-medium modifications are not accounted for so far. 



C. Medium modifications 

As medium modifications we implement, besides effects like Pauli-blocking, Fermi-motion of the nucleons and 
Coulomb forces, also different types of hadronic potentials for the particles and density dependent modifications of 
the decay widths. We emphasize that the imaginary parts of the potentials are in general not included since their 
effects are already accounted for by the collision term. Only in the case of the A resonance we follow a prescription 
which includes imaginary contributions of the optical potential as explained below. 

In position-space the nucleons are initialized according to Wood-Saxon-Distributions. For the momenta we utilize 
a local Thomas- Fermi- Approximation and express the local Fermi- momenta for the neutrons (pp) and the protons 
(Pf) as 

pJ'^'(f) = (37rV^(f))^ , 

where p"'P{r) denote the densities of the nucleons. The initial phase-space distribution of the nucleons is therefore 
given by 

fr:f,ir,p) = Ne{p'P^{r)^\p\). 

with N being a normalization factor. Due to the low particle multiplicities, the probability that a final state is 
Pauli-blocked PpB was approximated using the analytic ground-state expression 

p;^>-,p)^e(p7(r)-|pi) 

instead of calculating the Pauli-blocking probability based on ^""=1°°" as necessary in calculations of, e.g., heavy- ion- 
collisions. 



1. Baryon potentials 

The potentials are introduced as time- like components of vector potentials in the local rest-frame (LRF) ifTsh . For 
our purposes the most important mean field potentials are those acting on the nucleon, the A resonance and the tt 
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FIG. 5 The hadronic potential Aq of the A in the LRF as a function of momentum and density in units of MeV. 
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FIG. 6 The width Ftot of the A resonance a,t p — po in units of MeV as a function of momentum and invariant mass niA- The 
total width Ftot is the sum of the in-medium contributions and the Pauli-blocked free width F . 

meson. The nucleon potential is described by a momentum-dependent mean-field potential (for explicit details and 
parameters see (fish l Phenomenology tells us that the A potential has a depth of about —30 MeV at po 
Comparing to a momentum independent nucleon potential, which is approximately —50 MeV deep, the A potential 
is, therefore, taken to be 

Here we assume the same momentum dependence for the nucleon and the A potentials. The potential acting on the 
A is shown in fig. [S] as function of momentum and density in the local rest-frame. 

2. Modifications of the A width 

The A resonance is explicitly propagated in the model. Besides explicit iVA NA and NN A collisions 
and the vacuum decay channel A — > TriV, the model implements also an absorption probability for this resonance 
to account for A^A — > NN and NNA — > NNN processes This probability is based on the corresponding 
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contributions to the width of the A. derived in a microscopic model by Oset et al. jlj)- In fig. we show the total 
width 

Ttot ~ rmod,abs + (ATA + ^ vac 

of the A at normal nuclear matter density as implemented in the model. F denotes the vacuum decay channel, 
including Pauli blocking, and 

Tmod.abs — A >N N + ^NNl^ — >NNN 

is the absorptive contribution due to interactions with the nuclear medium as given by the model. 
3. Pion potential 

For the first time we have also included a realistic potential in the Hamiltonian of the BUU equation for the tt 
meson at very low energies. Below 80 MeV pion momentum a model by Nieves et al. l)38|) was used. In general, 
the pion spectral function shows at low energies three peak structures - the so called pion, nucleon-hole and A-hole 
branches. At momenta smaller than 150 MeV the pion is most probably found in the pion branch in which the 
dispersion relation does not differ much from the free one of the pion. Therefore we evaluate the optical potential at 
the vacuum position of the dispersion relation. 

In latest experiments at GSI t3S^ and MAMI |40l) new results on the pionic self energy at threshold were obtained. 
The results hint at a value of about A° (p = 0) ~ 25 MeV. The potential by Nieves et al. (|33) amounts to a somewhat 
lower value of about 18 MeV in symmetric nuclear matter at normal matter density. It is important to realize that 
such a repulsive potential limits the applicability of any semi-classical model to pion kinetic energies > 20 MeV in 
the central density region. 

Above 140 MeV a first-order A-hole result was used. In this energy regime the potential acting on the pion is defined 
by averaging over the spectral strength in the pion and the A-hole branches. To ensure a continuous derivative of the 
Hamiltonian, both models were matched in the region of p = 80 — 140 MeV. The potential for the pion is shown in 
fig. It is repulsive in the low-momentum regime, for higher momenta it turns attractive. The repulsion is caused 
by a only weakly momentum dependent S-wave term, while the attraction is due to a momentum dependent P-wave 
term, which vanishes at zero momentum and gains strength with increasing momentum. At constant density the 
potential is, therefore, continously decreasing with momentum in the range considered here. 
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FIG. 7 Hadronic potential Ao of the pion as function of nucleon density and pion momentum in symmetric nuclear matter in 
units of MeV. 
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III. RESULTS OF THE SIMULATIONS 

A. The mean free path of pions in nuclear matter 

Here, we discuss simple properties of the pion in the medium and compare to other theoretical calculations. In 
particular we investigate the total width 



= In 



rf(ln(iV.)) 
dt 



and mean free path A. These quantities are connected via the relation A = v/Ttot, where v denotes the velocity of 
the pions. In our model the loss of pions per time-step At is given by 



At - 2^ 

Jy —n,p 



E 

(Ar,M)e{(n.p),(n,n),(p,p)} 

where v denotes the pion velocity in the medium, ctjv - the nN scattering cross section and T^'^ is the absorptive 
three-body width of the pion and Ttot the total width. For the cross sections and the width we include the medium 
modifications discussed earlier. The integrals above are performed over the Fermi spheres of neutrons and protons. 

In practice, we obtain Ftot by performing a Monte-Carlo simulation with pions and nucleons initialized in a box with 
continuous boundary conditions within our BUU framework, including all the medium modifications and necessary 
collision rates. 

The momentum is calculated in the medium from the dispersion relation 



Etot^i = ^/pl+ml + A^{p) + Vc 

with the hadronic potential Aq and the Coulomb potential Vc ■ Finally, Ftot and A can be easily extracted by observing 
the time evolution of the pions. All results are shown as functions of the experimental observable 

^Vacuum _ Etotal ~ "T-ir 

in the classically allowed region £:Vacuum > A^{p) + Vc- 



1. Isospin dependence 



To begin with, we study the case in which we omit all potentials. In the GiBUU-simulation we use the same mass 
of TOjr = 138 MeV for all pions regardless of the charge. For both the positive and the negative pions the same data 
are used as input for the vacuum cross sections. Due to isospin symmetry the positive and negative charge states have 
the same properties in symmetric matter if electromagnetic forces are neglected. In fig.|Slwe observe that the positive 
and neutral pions do actually have the same velocities, but not the same width in nuclear matter. The difference is 
quite small except for very low energies. The reason for this can be found in the elementary absorption and scattering 
processes. As mentioned before, we have used vacuum data to pin down the elementary reaction rates. Therefore the 
absorptive width of the charged pions is fixed by the elementary reactions 

p p — > p n TT^ 
p n ^ p p Tr~ 

while the absorptive width of the tt" is determined by 

P P ^ p p TT^ 

p n p n i:'^ 

through detailed balance. Fig. 0] shows that the three-body absorption rate for the neutral pions at low energies is 
higher than for the charged ones. This explains the larger width of the 7r° near threshold. 
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FIG. 8 BUU results without electro-magnetic forces inside symmetric nuclear matter at po ~ 0.168/m~ . The results are 
shown with (dashed) and without (solid line) hadronic pion potential. The upper panels show the velocity of the pious in 
nuclear matter, the lower ones visualize the full width of the pion. The lower curves in the lower panel denoted "only 3-body" 
do not include two-body processes or resonance production. 



2. Collisional width in the medium 



Using the neutral pion as an example, we first discuss the lower right panel of fig. |H1 At low energies there 
is practically no difference in the width between the simulation with (dashed lines) and without (solid lines) pion 
potential. Due to the absence of resonance contributions in this energy regime, the width of the pion is dominated 
by its absorption via the NNt: — > NN background. At higher energies a small difference in the two curves can be 
observed. Here the decay width is somewhat smaller if the pion potential is included. 

This difference can be understood by investigating the resonance production process Ntt ^ A in detail l|20t) . At 
fixed pion energy, the attractive pion potential leads to a higher momentum of the pion and, therefore, also to a 
higher momentum of the produced A, resulting in a decrease of its mass. In the considered energy regime this mass 
is smaller than the pole mass. As a consequence, the production of the A resonance is less probable and the width of 
the pion decreases. 



3. Velocity of the tt mesons in medium. 

The velocity of the pions is given by Hamilton's equation 

dn _ dH__ d(^^p^+m-^ + Al + Vc) ^ ^ dA^ ^^^^ 

dt dpi dpi ^ ^2 dpi ' 

The hadronic potential is repulsive for low values of Et^t and of the same order as -Bkin- Therefore, in the first term of 
equation 110|1 the pion momentum is considerably smaller and, as shown in the upper panels of fig. |H1 the velocity is 
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FIG. 9 Mean free path of a neutral pion in symmetric nuclear matter at po — 0.168/m~'^. We show the result without 
(solid line) and with (dashed line) hadronic potential for the pion; all other standard medium modifications are included. For 
comparison, the result of Mehrem et al. (|4lh is also shown (crosses linked by a dotted interpolating line), which will be discussed 
later. 
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FIG. 10 Influence of the electro-magnetic forces on the mean free path of the charged pions (p — 0.168fm~^ , pp = pn = f ). 
On the left panel we show results of our model without Coulomb potential - in this case the charged pions have identical 
mean free paths. The curves presented with dashed dotted lines are charge-averaged results of P. Hecking di^) . which we will 
discuss later. The panel on the right shows the mean free path of the charged mesons, including Coulomb effects. We explicitly 
investigate the result with and without hadronic potential and Coulomb potential for the pion, all other standard medium 
modifications are included. 



strongly modified at low values of Ekin- The second term dA^ /dp is always negative. It is large in absolute magnitude 
for low energies and small for energies greater than 80 MeV and leads to an overall reduction of the velocities. 



4. Results for the mean free path. 

After having extracted the width of the pion in the nuclear matter restframe directly from our numerical simulation 
the mean free path is obtained by A = v/T. Therefore we now consider the two effects discussed earlier in order to 
understand the changes in the mean free path: modification of the decay width F and modification of the velocity v 
due to the potentials. The effect of including the hadronic potential becomes visible in fig. O especially at very low 
energies. The mean free path drops rapidly at low energies compared to the simulation without hadronic potential. 
This sharp decrease of the mean free path as a function of asymptotic kinetic energy stems from the velocity decrease 
at low energies when including the repulsive hadronic potential. At larger values of the kinetic energy the effect of 
width and velocity just compensate each other. 

For charged pions, electromagnetic forces play an important role, especially at low energies. Since in nuclear matter 
calculations we cannot introduce any Coulomb potential based on the actual density, we assume a reasonable constant 
Coulomb potential of Vc = ±20 MeV depending on the charge of the pion. The results for the mean free path are 
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FIG. 11 The mean free path of a tv^ investigated at different densities of the nuclear matter. The hadronic potential for the 
pion is not included. 



shown in fig. ^1 The most important effect of the hadronic potential is the large decrease of the mean free path at 
very low energies, qualitatively similar to that already discussed for the tt*^ above. 

We have also studied the density dependence of the mean free path, choosing the tt" as a showcase. The results are 
shown in fig. 1111 It is important to note that the density dependence is highly nonlinear, contrary to the low-density 
limit. This non-linearity is generated by the NNt: — > NN process, which goes to first order quadratically with density, 
and by the implicit density dependence in the medium modifications. 



5. The vacuum approximation. 

We now compare our results to the eikonal vacuum approximation l)35l: ksTl .commonlv found in the literature. In 
this approximation the mean free path is given by 

2 is[p) p avac 

denotes the imaginary part of the momentum. Using the total cross sections shown in fig. El we can now 
evaluate the mean free path. We assume cr^-n^x = '^tt+p^x a-nd cr^+„^x ~ (^-^-p^x because there exist no data 
for TT scattering off neutrons. 

In the case of symmetric matter, we calculate <Jvac for the tt"^ and 7r~ by averaging over the neutron and proton 
contribution to the cross section: 

^Vac = 2 ('^^^"^^ + cr^±p-^x) 2 {^-^^p-^x + <^7T-p^x) ■ 

For charged pions the mean free path in this approximation is shown in fig. 1121 Comparing it to a full BUU calculation, 
which considers the Pauli blocking of the final states, as well as the Fermi motion of the initial states and three-body 
interactions, we see a dramatic discrepancy at small energies. 

For completeness we have also investigated the influence of Fermi-motion and calculated the mean free path by 
integrating over the interactions in the Fermi sea: 

^ = ^- TT- ■ (12) 

PF 

The result can also be seen in fig. 1121 The modifications of the mean free path due to the Fermi motion in the vacuum 
approximation are very small. 

The vacuum approximation does not consider any kind of many-body interactions at low energy. Only the two body 
■kN interaction is allowed. In our model, detailed balance dictated us to introduce a NNtt NN process to account 
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FIG. 12 Mean free path of the charged pions in the naive vacuum approximation. The inclusion of Fermi motion does not 
have a sizeable impact. A BUU calculation without Coulomb forces and without hadronic potential for the pion is shown for 
comparison. 



for the NN NNtt channel. In fig.|Hlwc sec the large contribution of this process to the width at low energies. Such 
a three body NNtt interaction is not included in the vacuum approximation. At higher energies, especially above 
120 MeV the vacuum approximation coincides with our results, since there the dominant process is the two body 
process Ntt — > A. 

We conclude that the naive vacuum approximation is qualitatively and quantitatively not reliable in the energy 
regime of Skin ^ 70 MeV where multi-body collisions, potential effects and Pauli-blocking are important. 



6. Comparison to optical model results. 

In lj43l Hecking published calculations of the pion mean-free-path based on two different types of phenomenological 
optical potentials Vopt- The first parameter set for an optical potential is taken out of Strieker's analysis 13), the 
second one is calculated by Chai and Riska l|45t) . Starting from equation Hll|) one can calculate the imaginary part 
of the momentum using the dispersion relation p'^ = — {m? + 2E Vopt{p,E)) approximating 11 = 2E Vopt by 
ignoring terms of the order O (V^p^/m^). Instead of solving the real and imaginary parts of the dispersion relation, 
Hecking approximates the real part of the dispersion relation by the vacuum solution [p^] = E^ — rn? and uses this 

approximation in the equation for the imaginary part 5 [p^] = —3 2E Voptip^E = \/p'^ -\- m?) ■ This last equation 

now defines also 3 [p], which can be used in equation l|ll|l to obtain the mean free path. To minimize Coulomb effects, 
he defines an averaged mean free path 

^avg = 2 + ^TT- ) • 

In order to compare to his results we have performed a calculation without Coulomb forces. In fig. Hecking's 
results are compared to the full BUU calculations. As a most prominent feature one observes that our results are much 
lower than Hecking's over the full energy range. The work by Mehrem et al. ()4ll) explains this discrepancy. Solving 
the full dispersion relation with an optical potential calculated by J. A. Carr ^, it reports a mean free path of the pion 
which is smaller than Hecking's result. In this work the large difference is attributed to Hecking's approximations 
described above. 

As a benchmark for our model, the result of Mehrem et al. (|4ll) for the mean free path is shown in fig. |51 in 
comparison to our result obtained with the BUU simulation ^. Including the hadronic potential in our model we 
notice that the mean free path decreases considerably at low pion kinetic energies. There the hadronic potential 



^ The parameters are given in i4ll) 

^ The real part of Vopt, which is used in Mehrem's work, is very strong at low energies with Vopt ~ 48 MeV at E^f^ = 0. A direct 
comparison by including this real part in our semi-classical model is, therefore, not possible. 
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becomes repulsive with V > E^in and the semi-classical model breaks down, as mentioned earlier. In contrast to this, 
quantum mechanics allows for tunneling, i. e. propagation, into such classically forbidden regions. 

A proper discussion of the mean free path is obviously important in the analysis of experiments with final state 
pions being produced inside the nuclear medium. When working in a quantum mechanical framework solving the full 
dispersion relation, as done e.g. by Mehrem et al., is very important. Ad hoc assumptions have a large effect on the 
pion mean free path and this has a dramatic effect on any final state analysis. 

Since the mean free path is not directly observable, it is ultimately an open question whether transport gives a 
reasonable mean free path for the pion. This can only be answered by experiment. A test for our model assumptions 
will, therefore, be absorption cross sections which we address in the next paragraph. 



B. The scattering of pions off nuclei 



Low-energy pion scattering experiments have been studied extensively with elementar y ta rgets ( e.^ 
however, there exist only a few data points for pions scattering off complex nuclei (|47l : l48i l49t ISOt ISU 1521 : \5c 

The total cross section (f5ll : IH^I) includes the coherent contribution which would have to be substracted, before 
comparing them to our results (for a discussion on this issue see e.g. (ji^)). We study here solely reaction, charge 
exchange (CX) and absorption cross sections. In the considered energy regime, reaction and CX cross section mea- 
surements are rare (^3 l4^ . In fig. El we present our results for reaction and CX cross sections for ^^C and ^^^Bi. 
The overall agreement with the few existing data points is satisfactory. Fortunately, the experimental situation for 
absorption cross sections |4^l49H50|) is more promising. Therefore, this observable can be used to evaluate the quality 
of our model. 
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FIG. 13 Reaction ("re ac" ; solid lines) and charge exchange ("CX"; das hed lines) cross sections for ^^C and '^°^Bi. The data 
points are taken from Ji^ (open circles: acx, open squares: (Treac) and (full quares: (Jroac). 



In a quantum mechanical approach one cannot calculate absorption cross sections in the general case. Only at very 
low energies, where the quasielastic contribution to the total reaction cross section is negligible, such an absorption 
cross section can be obtained(e.g. (fssl) ). One thus relies on eikonal approximations to split the reaction cross section 
into a quasi-elastic and absorptive part. On the contrary, in the BUU simulation the calculation of an absorption 
cross section is straightforward due to the semi-classical treatment. 

In figures 1141 and 1151 we show calculations for different nuclei and different bombarding energies of the pions. 
Comparing the curves obtained without any potential (dotted) to those with the Coulomb potential included (dashed), 
we see that the Coulomb potential alone has only a small influence on light nuclei (fig. I14|l . but is very important at 
low energies for heavy nuclei (fig. I15|l . Its long range leads to a sizeable deformation of the trajectories already long 
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before the pions reach the nucleus. Therefore the negative pions can interact with the nucleus even if they have large 
impact parameter while the positive pions are deflected. In presence of the Coulomb potential we see a reduction 
of the cross section for the 7r+ and a large increase of the cross section for the tt^. This agrees to the findings of 
Nieves et al. l)54j) . who pointed out the relevance of the Coulomb potential in their quantum mechanical calculation 
of absorption and reaction cross sections. 

When one includes the hadronic potential for the pion, another overall effect sets in. Once the pion enters the 
nucleus it is influenced by the short-range hadronic potential, which amounts to —40 MeV at high momenta and to 
+20 MeV at low momenta, as well as the Coulomb potential which amounts to roughly ±10 MeV in a medium size 
nucleus, and to roughly ±20 MeV in the case of Pb. At very low energies the two potentials nearly compensate for 
the negative pion, while they add up to a strongly repulsive potential in the case of a 7r+ . 

In fig. El we see a fair agreement on C and Al with the data, when the hadronic potential is included. Addressing 
heavier nuclei, fig. 1151 shows the absorption cross sections for Cu, Sn and Au. The overall agreement to data is very 
good if the hadronic potential is included, especially for Sn and Au. The curves which include a hadronic potential 
for the pion show a prominent kink structure at roughly 30 MeV (20 MeV) for the positive (negative) pion. This kink 
is caused by the repulsive character of the pion potential at low energies. On one hand, the mean free path decreases 
rapidly at very low energies (compare fig.EJ - this causes the absorption cross section to rise. On the other hand, the 
repulsive potential pushes the pions out of the nucleus or even reflects them. Below the kink the repulsive feature is 
more prominent; above the kink the decreasing mean free path is more important. 

As a conclusion, we find that it is critical to include Coulomb corrections. On top, the absorption cross sections 
are sensitive to the hadronic potential of the pion, i.e. to the real part of the self energy in the medium. As we have 
already seen in fig. ^1 the mean free path is quite insensitive to the hadronic potential except at very low energies. 
We thus conclude that the modification of the trajectories of the pion is the main effect of the hadronic potential. 
In its repulsive regime the hadronic potential pushes the pion outwards and the overall path of the pion inside the 
nucleus becomes shorter. Therefore the probability of absorption is decreased. The attractive behavior at larger 
energies causes the opposite effect. 

The overall agreement to the data is satisfactory in spite of some discrepancies, especially in Al and Cu. Considering 
the fact that the pions have very large wave lengths at such low energies, the success of the semi-classical BUU model 
is quite astonishing. Due to the large wave length one expects also many-body correlations and quantum interference 
effects to be important. Many-body effects arc partially included via the mean fields acting on pions and baryons 
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FIG. 14 Absorption on light nuclei depending on the choice of potentials for the pion. The data points are taken from 
Ji^: ligtlsol) . In all plots the standard hadronic potentials for the nucleon and A resonance are used. 
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FIG. 15 Absorption on heavier nuclei depending on the isospin of potentials for the pion. Data points are taken from l|4flll . 
In all plots the standard hadronic potentials for the nucleon and A resonance are used. We explicitly vary the inclusion of 
Coulomb potential and hadronic potential for the pion. 



and the modification of the A width. Besides this we included only 1 <-^- 2, 2 ^ 2 and 2^3 body processes in the 
collision term. We take the success as an evidence that no higher order correlations than the latter ones are necessary 
to describe pion absorption. 



C. Photon-induced tttv production in nuclei 

1. Overview 

We now discuss the consequences of the pion mean free path on recent double pion production data off nuclei^hich 
have been interoreted as a source of information about the spectral density in the a channel in nuclear matter (|9t l57f) . 

Already in l|13^ we have presented results on the double-pion photoproduction off nuclei using the transport ap- 
proach. There the treatment of the pion final state differed from the treatment presented in the previous sections. 
In order to compare to a work by Roca et al. l)57|) . we implemented there the same absorption probability for the 
pion as in their work. Therefore we did not propagate the resonances explicitly, but rather used experimental data 
to describe the quasi-elastic Ntt — > Ntt scattering. We also used contributions to the imaginary part of the optical 
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FIG. 16 Two pion invariant mass distributions for tt^tt" and tt^tt" photoproduction ofT tlie proton and neutron obtained from 
the model of iSStlSGt) . We show results on two different energy bins, E-y = 0.4 — 0.5 GeV (upper curve) and Ej = 0.4 — 0.46 GeV 
(lower curve). Data for Ej — 0.4 — 0.46 GeV are taken from 



potential of l)38() to model pion absorption. The aim there was to indicate the relevant role of conventional final-state 
interaction (FSI) effects in two-pion photoproduction in nuclei. 

The authors of l)57|) achieved quite an impressive agreement with the experimental data, when studying double pion 
production in a many-body approach which allowed for pion-pion correlations. However their treatment of the pion 
FSI was based upon a purely absorptive Glauber damping-factor calculated along straight line trajectories. 

In the previous section we have shown that our transport approach is still reliable in the regime of long pionic wave 
lengths. Since there is no fully quantum description for the incoherent 2tt reaction available, we utilize the semi- 
classical final state description to simulate the TTTr-photo production experiment performed by the TAPS collaboration 
ifol: ITol: lilt IT^ . Originally this experiment was intended to concentrate on the lower 400 MeV < < 460 MeV bin, 
focusing on the threshold behavior of the tttt channel. Lately the experiment has evaluated a larger energy bin for the 
photon energy in order to enhance statistics. Therefore, also results for the larger bin with 400 MeV < < 500 MeV 
are presented here. 
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FIG. 17 The cross section for vr tt and vr tt photoproduction off the proton at E-y — 0.45 GeV. The mass of the 

pion-pair is denoted m^^ whereas p denotes the absolute momentum of one of the pions. 



2. Model for the elementary reaction 

For the elementary two pion production process on the nucleon we exploit the model of 1)551 : [^^. This model 
provides a reliable input for the mass distributions of the pions in the elementary process. The total production 
cross sections are chosen according to the model in those channels where the agreement to data is satisfactory 
ilP ~^ n7r"'"7r°,7n riTr+Tr"). In the remaining ones the model predicts lower total cross sections than those observed 
in experiment at the ene rgy of interest for us. Hence, we use directly the data measured by the TAPS and the 
DAPHNE collaborations ifssl : Isol : l60l : Isil l63l to normalize the calculated cross sections, while we take the decay 
mass distributions from theory. The implemented elementary invariant-mass distributions are presented in fig. 1161 
for both a neutron and a proton target. Presently, only data for invariant mass distributions on a proton target are 
available. The solid curves in fig. El represent the results of the elementary model (j^ IHfil) for tttt production for 
E~f = 400 — 460 MeV compared to the experimental data. The model describes the shape of the data quite well; 
for the semi-charged channel the cross section is somewhat underestimated. The dashed curves in fig. 1161 show the 
predictions for = 400 - 500 McV. 

To illustrate the momentum distribution of the pions, we concentrate on a jp reaction with a fixed photon energy 
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FIG. 18 Two pion invariant mass distributions for tt^tt" and tt^tt'' photoproduction ofT ^^C, ""'Ca and ^°®Pb for = 0.4 - 
0.46 GeV. The uncharged tt'^tt" channel is plotted on the left, the charged channel tt^tt" is presented on the right hand side. We 
show results without final state interactions (dashed dotted) and results with FSI and respectively with (dashed) and without 
(solid) hadronic potential for the pion. 



of 450 MeV. This distribution is shown in fig. El function of both mass of the pion-pair and momentum of one 
of the pions. Most of the pions are produced with momenta around 120 MeV, corresponding to iJkin ~ 44 MeV, i.e. 
just at the lower hmit of the apphcabihty of our semi-classical model at pQ. Note that the momenta of the produced 
pions range between — 250 MeV. In nuclei, the FSI is therefore dominated by the A resonance, S'-wave scattering 
and the NNtt NN process. 

3. Results and discussion 

Our results for tttt photoproduction off C, Ca and Pb nuclei are presented in fig.^jand fig. ll9[t he two shown energy 
bins are the ones which are presently considered in the ongoing analysis by the TAPS group ijgt ITll IT^ . As already 
discussed in ((l3|) . we observe that absorption, elastic scattering and charge exchange processes cause a considerable 
change of the spectra with the peak of the mass distribution moving to lower masses due to rescattering. 

The hadronic potential for the pion causes only a minor effect, so the observed reaction is rather insensitive to this 
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FIG. 19 Two pion invariant mass distributions for tt^tt" and tt^tt'' photoproduction ofT ^^C, ""'Ca and ^°®Pb for = 0.4 - 
0.5 GeV. The uncharged channel is plotted on the left, the charged channel tt^tt" is presented on the right hand side. We 
show results without final state interactions (dashed dotted) and results with FSI and respectively with (dashed) and without 
(solid) hadronic potential for the pion. 



aspect of the model. This can be understood by analysing the production points of those pion pairs which are not 
absorbed and which are, finally, observed. In fig.lJOlthe cross section for tt'^tt'^ production off Pb at 500 MeV is shown 
as function of the production point R of the pair. Without FSI the distribution da^o^o / dR is proportional to p{R)R^; 
including FSI the distribution is shifted and centered around 6.8 fm which corresponds to roughly p — 0.076 fm""^. 
Hence, the potential has a minor effect on the observed pions due to the low density at the initial production point. 

Since the mean free path increases at lower momentum, a cut that keeps only low-momentum pions could reduce 
the sensitivity of the whole production process to the 7r-nucleus FSI. On the other hand, this would complicate 
the theoretical treatment of the process considerably, since quantum mechanical effects would become even more 
important. 

First experimental results for this process have been presented in (Q). Recently, those data have been reanalyzed 
by the TAPS collaboration. The preliminary results of this analysis can be found in |Tlt IT^ . and a comparison 
of the former results to this new analysis is given in l)lO() . The old analysis is being re-evaluated, especially in the 
semi-charged channel. Therefore, we do not present these data sets here. At the moment, new sets of data, taken 
with a An setup with Crystal Ball and TAPS at MAMI, are being analyzed. With these a meaningful comparison 
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FIG. 20 Calculation without hadronic potential for the pion for Ej = 500 MeV off Pb. The cross section for tt^tt" production 
is shown with FSI (solid steps) and without FSI (dashed steps) as function of the original production point R of the pion pair. 
For illustration we show the density p of the Pb nucleus (dashed dotted curve). 
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FIG. 21 Illustration of the uncertainty due to the elementary tttt production process. We show standard calculations without 
hadronic potential for the pion for = 400 — 500 MeV. Only the cross sections used as input for 7A'^ — > Nmr are varied 
between experimental data (solid lines), experimental upper (long dashed lines) and experimental lower bounds (short dashed 
Unes) of all cross sections. 



with experiment will be possible. 

The main theoretical uncertainties in our calculation stem from two different sources - pion reaction mechanism and 
initial production cross section. The pion FSI reaction mechanism incorporated in the model and its uncertainties 
have been discussed in the preceding sections. The influence of the uncertainty in the elementary cross sections is 
illustrated in fig. 1211 The three curves shown in each plot use as input either the given data points or, simultaneously 
for all data sets, the lower and upper bounds. One has to note, that jn — s- p7r~7r'^, jp — > p7r"^7r~ and especially 
jn —> TiTT^Tr'^ l|60() are sources of uncertainty in the elementary cross sections. We emphasize, that this scenario of 
simultaneous lower and upper bounds instead of a computational intensive x-squared-analysis is only meant to give 
an estimate for the error. The statistical error due to the Monte-Carlo method was estimated to be about 2% by 
comparing parallel runs, being therefore negligible. 

To summarize this section we emphasize that final state interactions of the pions are strong and tend to shift the 
maximum of the tttt mass distribution in all channels towards lower masses. This effect considerably complicates 
drawing a link between the experimental data and a possible softening of the in-medium / = channel. Any 
theory aiming to describe the observed effect on the basis of a partial chiral symmetry restoration or an in-medium 
modification of the tttt production process must take the final state effects into account. 

IV. SUMMARY 

In this work we have presented the treatment of low-energy pions within the GiBUU transport model. In particular, 
we have studied the mean free path and compared to optical model results. Studying the long wave-length region. 
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we have pointed out the problem with pion energies below 20 MeV due to the semi-classical treatment in our model. 
As a benchmark, we have investigated pion absorption cross sections with very good agreement to experimental data 
from light to heavy nuclei. At low energies, we have found a strong dependence of this observable on the Coulomb 
potential and the pion hadronic potential. Finally, we have shown predictions for jA — > ttttA in the regime of 
400 MeV < < 500 MeV and also new calculations for 400 MeV < < 460 MeV. We emphasize, that the link 
from TTTT photoproduction experiments to imderlying effects like a possible softening of the in-medium tttt interaction 
is considerably complicated by the necessary presence of pion final state interactions which also shifts the tttt invariant 
mass distributions towards lower masses. 
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